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ABSTRACT 


This work concerns numerical studies of PVT air collector operating in 
dynamic mode. We presented the temperature of the solar cells and of the air 
at the outlet of the PVT collector with and without reflectors. The obtained 
simulation results have been validated by comparating with the experimental 
results available in literature, where good agreement was been noted. In 
addition, we determine the optimal position of the reflectors by numerical 
calculation, in order to obtain the maximal concentration of the solar radiation 
intensity on PV/T collector. The thermal and electrical efficiency of PV/T 
collector without reflectors and with reflectors in optimal position have been 
calculated for composite climate at Tetouan (-5°N, 35.5°E). So, the effect of 
adding solar reflectors on the electrical and thermal performance is showing. 
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Solar energy 
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1. INTRODUCTION 

Solar energy, a type of renewable energy, has a positive effect on the environment; contributes to 
minimizing CO? emissions into the atmosphere. Currently, it is operated by a variety of applications, including 
space heating, water heating, drying of agri-food products, solar refrigeration, and electrical energy produced 
by systems photovoltaic. 

The integrated of the solar thermal and solar photovoltaic systems have evolved considerably, the two 
technologies provide the opportunity to increase use of solar energy and enhance efficiency, these systems 
PV/T collectors deliver both electrical energy and thermal at the same time, also provide more useful energy 
of collector than separate photovoltaic and thermal systems. 

Some of the earliest research papers on PV/T systems date back to the 1970s, Kern Jr, E.C., and 
Russel, M.C [1] and Florschuetz L.W., [2] have recently received renewed interest and expected to become 
increasingly interesting in the future, because the combination of the two technologies offers the opportunity 
to increase efficiency and expand the use of solar energy. 

A bibliographic synthesis of PVT solar collectors was carried out, which enabled us to note that 
research on solar collectors have been the focus on improving performance and innovation, including the work 
of A. Tiwari, M.S. Sodha, A. Chandra, and J.C. Joshi [3] which investigated the performance of the 
photovoltaic (PV) module integrated with air duct. They validated the theoretical and experimental result for 
composite climate of India and they concluded that an overall thermal efficiency of PVT system is significantly 
increased due to utilization of the thermal energy from PV module. A.S. Joshi, A. Tiwari, G.N. Tiwari, I. 


Journal homepage: http://ijpeds.iaescore.com 


1046 O ISSN: 2088-8694 


Dincer, and B.V. Reddy [4] have studied the performance of unglazed hybrid photovoltaic thermal (PV/T) 
glass-to-glass system for composite climate of New Delhi and the comparison with glass-to-tedlar (PV/T) 
system. 

S. Kumar, and A. Tiwari [5] Investigated the concept of a novel self-contained integrated 
photovoltaic/thermal (PV/T) hybrid solar still, which was developed and verified for a composite climate at 
[.1.T. New Delhi. This PV system has been employed to generate electricity to run the pump (60W and 18V); 
the innovative structure of the hybrid-active solar still allows working in any remote location due to its 
autonomous nature. 

Sandnes, B., and Rekstad, J. also [6] tested the PV/T unit experimentally to determine thermal and 
photovoltaic performance, and used a polymer solar heat collector that has been combined with 
monocrystalline silicon photovoltaic cells in a production unit of hybrid energy that simultaneously produced 
heat and electricity at low temperatures. As well, Y. Tripanagnostopouls [7] studied hybrid PV/TH solar 
collectors with either air or water as a heat transfer fluid, which can be integrated into buildings. The objective 
of their work was to reduce the operating temperature of the PV modules, to increase the production of 
preheated air and to reduce thermal losses through the insulation. 

In addition to that, Ahmad Fudholi, et. al. [8] performed analytical studies on energy and exergy 
efficiency of the PVT air collector with monofacial and bifacial solar cells, the latter is made up of two active 
surfaces that can capture more of the solar radiation with their front and back surfaces, therefore more 
performance compared to conventional monofacial PV. P. Dupeyrat, C. Menezo, and S. Fortuinc [9] compared 
numerically and experimentally a photovoltaic-thermal system with a conventional solar system, they indicated 
that (PVT) is the more advantageous system, with respect to energy saving and primary energy consumption. 
N.S.B. Rukman, and A. Fudholi [10] Presented the energy and energy efficiency of water-based PVT systems 
and summarized previous studies on the development of this system, also explained water PVT system models. 

R. Kumar, and M. Arosen [11] Have critically review photovoltaic-thermal solar collectors for air 
heating and recent advances, their discussions Included of the status, development and applications of various 
PV and solar thermal technologies, and major developments in PV/T collectors and PV/T air heater systems. 
Their objective is to increase understanding of solar PV/T collectors and heaters and help increase their 
adoption. Y. Yu, E. Long, Xi Chen, and Hongxing Yang [12] have installed Unglazed Photovoltaic Thermal 
Collector in Chengdu, Sichuan Basin were input the measured ambient data and operation data to two TRNSYS 
models and were compared the results, by calibrating the models with experimental data, they reduced the 
deviation between the simulated and the measured values. In addition, the excellent performance of the roll- 
bond PVT absorber was also validated by the steady-state test results. 

As for studies related to flat reflectors, LJT. Kostic, TM. Pavlovic, and ZT. Pavlovic [13] reported on 
the results of the influence of reflectance of flat-plate solar reflectors, made of aluminium sheet and aluminium 
foil, on the energy efficiency of PV/T collectors. They presented both experimental and analytical results of 
the optimal position of aluminum sheet made flat plate solar reflectors during the daytime over the whole year. 
Similarly, H. P. Garg, R.K. Agarwal, and Ashok Kumar Bhargava [14] studied a PVT-air hybrid system that 
included a planar booster to increase the intensity of sunlight on the photovoltaic module. They calculated the 
minimum solar cell area required to operate the pump at a specified rate without and with the boosters. 

In this paper, the principal objective is to study the influence of the two planar reflectors on the thermal 
and electrical performance of a PVT hybrid collector located in southern Morocco. We carry out numerical 
simulations to determine the optimal inclination angles of the reflectors, and calculate the temperature of the 
air and the cells without and with reflectors at this site. 


2. HYBRID PVT AIR COLLECTOR 
2.1. Description of the studied system 

Photovoltaic thermal (PV/T) collectors are devices that simultaneously convert solar radiation into 
electricity and heat. the simplified diagram of the photovoltaic thermal hybrid collector PV/T collector shown 
in Figure 1, It consists of a glass cover, monocristalln’s photovoltaic module of 1.29m length and 0.33m width, 
aluminium absorber placed at the lower part of the module, this absorber is enclosure containing the coolant 
and insulation to reduce wool heat loss. 

The system shown in Figure 2 consists of three different components: the PVT collector, the upper 
reflector and the lower reflector. The collector is tilted at 35° and faces south. Both reflectors are tilted to the 
horizontal plane by al and to the vertical plane by a2. 
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Figure 1. Schematic of PV/T model 
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Figure 2. Photovoltaic thermal collector with two reflectors 
2.2. Global thermal expressions for the PV/T solar collector 
The heat balances for each node, in a section slice, are written according to the general shape: 


c; dT; 
As dt 





= P, + $ (hci + hrij + hij)(T, —T,) (1) 


The energy balance for the different parts of the PVT solar collector according to [15]. 
For the glass cover is: 


M,C 
et (T; -Ty) =P, +hy (Tọ Ty) + Bigot (Tax -Ty ) + Beg -cet (Toen -Tg ) +8 g-cet (Ten -Ty) 2) 


for the photovoltaic cells (PV) is: 
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where: 
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for the absorber plate 1s: 


ppp 
AsAt 


a oe | 


T, -Tp =P, +becen-p (Tee -Tp )+ hp-t (Te -T )+ Bopp (Tos - Tp] (5) 
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alr Gat -To | =P,, +h, ¢(T;-Ty,)+h 


AsAt P (T, = Th F Heppin (Ta g T ) at Dhan (Ta =L ) (7) 
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for the insulation 1s: 


mah (Ta 7 Ta ) z Dba (T =Ta ) $ Rejoin (Tp = dig )+ Din-s (T, -Tn ) + Bring (T, -Tn ) 8) 


2.2.1. Convection heat transfer coefficient 


The convection exchanges that occur in the system under study are: 
a. The convection between the glass and the atmosphere, calculated by Mc Adams [16]: 


h, =5.67+3.86V,, (9) 


b. The natural convection of air between the glass and the PV model 


À, Nu 
—cell 7 
i e 


h 
g—cell (10) 


c. The convection between the fluid and the absorber plate 


r 
hp- = hp-t = Nuf a (11) 


h 


Nuf Is the Nusselt number is calculated according to the flow regime; however, the Nusselt groupings 
corresponding to these regimes are expressed by the equations in [17]. We took the condution and the 
coefficient of radiation heat transfer from [18], [19]. 

To solve the equations of this model, the iterative numerical method used is based on the method of 
Gauss seidel [20]. Indeed, a program was written in Fortran 90 in order to find the values of the temperatures 
Tg, Tcell, Tp, Tins and Tf of the air collector, as well as the value of the thermal efficiency. 


2.2.2. Thermal efficiency 

The thermal efficiency of the PV/T collector can be expressed as a ratio of the amount of heat extracted 
by the fluid used, denoted Qu (W/m2), to the amount of solar radiation incident on the glazing. The thermal 
efficiency is determined from [21]: 


mC, (T; -T,) 


D — sea (12) 
A, Gio 


2.3. Electrical efficiency 
The electrical efficiency Npy of the crystalline silicon (c-Si) solar cells/modules presented by [22]: 


Nov = Net [1+ BCTeen — 298)] (13) 


Where 77,¢,and p respectively represent the coefficient of photovoltaic conversion efficiency and the 
coefficient of photovoltaic conversion efficiency at reference temperature (298 K). 


2.4. The reflectors equations 

The theoretical equations that allow us to carry out the algorithm to calculate the global radiation 
received by inclined PVT collector, are the formula of the solar altitude and azimuth angles, also we use the 
program FORTRAN as a tool of test and simulation to determine the trajectory of the sun immediately. 
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The instantaneous solar altitude angle h is given by [23]: 
sinh=singsind+cos@cosdcos@ (14) 


where: 
© Is the latitude of the location 


0 Is the declination angle given by. 
0=23.54°sin[0.980°(N+287) | (15) 


w Is the hour angle, which defined as: 
o=15°(TSV-12) (16) 
The extraterrestrial beam radiation reaches the ground is: 
AMO-078 


i=i,y, *0.7 (17) 
With 


-1 
AM=] cos0,+0.50572(96.07995+0, D (18) 


The approximate intensity of solar radiation on the ground at a given location can be measured by 
calculating the extraterrestrial intensity ig on a surface that is perpendicular to the sun on a given day of the 
year, according to N. Kumar, G.Vishwanath, A.Gupta [24]: 

loh =o “cos0, (19) 


toiset oonasreos{ 2D) (20) 
365.25 


Where i,.= 1367 W/m? is the solar constant, and 0, is the zenith angle of the sun, correlated to the altitude 
angle as (21). 


6, =90-h (21) 


In order to calculate the solar radiation G on the inclined PVT collector with B, we have used the 
following expressions [25]: 


G=Gdirt+Gdif+Gert+Gr1+Gr2 (22) 
Where: 
Gdir=Gin*sin(B+h) (23) 
1 
Ga OP (24) 
2 
1-cos 
Ger=pe*Gh* mat (25) 
Grl=pal*Gin*siny*sin(h-a1) = y=B+2*a1-h (26) 
Gr2=pal*Gin*sint*sin(h+o2) => t=ht+2*a2-B (27) 
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3. RESULTS AND DISCUSSION 

The total radiation falling on PVT collector tilted at a fixed orientation with and without reflectors, 
and the ambient air temperature are plotted in Fig.3. The collector is oriented 35.5° to the south, while the angle 
of the reflectors can changed al the angle of the lower reflector is changed from 7° to 82°, while the angle of 
upper reflector is changed from 10° to 35°. So, we optimize the angles of the reflectors to reach a higher 


intensity of solar radiation falling on the PVT surface, from the curves of Figure 4, we conclude that the optimal 
value of al is 37° and a2 is 10° during May 2 at 12:00 noon. 
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Figure 3. Hourly solar radiation on PV/T collector Figure 4. A change of total solar radiation on PV/T 
with and without reflectors collector depending on reflectors position 


In Figure 5, we have presented the air outlet temperature given by Touafik [26] and our model without 
reflectors in June 15. We observed that the result of our simulation is in good agreement with the experimental 
results. 

In this comparaison, we kept the same properties of the system used in the experiment, except the 
absorber plate equipped with rectangular fins. Although not used the fins. The air outlet temperature of our 
model is almost equal to the temperature of [26], so we have found that our model is efficient can be used to 
improve the performance of the PV/T collector. 

Figure 6 show the air temperature at the outlet of the PVT collector for 4 may without and with 
reflectors in optimal position, we notice that the latter is higher than the outlet temperature of the air in the 
collector without reflectors. Because these temperature values are related to the intensity of solar radiation. 
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Figure 5. Validation of the mathematical model of Figure 6. Effect of reflectors on the temperature 
Hourly variations of the outlet air Temperatures of Outlet of air of PVT collector 


the PV/T collector 


Figure 7 illustrates the comparaison of the temperature of the solar cell without reflectors and with 


reflectors of our numerical program with the temperature without reflectors available in the experimental study 
of [27] (4 May 2016). According to the graphs, there is a concordance between their experimental results and 
our theoretical results. The temperature of the solar cell varies from a minimum value of 320 K at 8h00 to a 
maximum of 340 K at 13h00, with a root mean square between the experimental values and those of our 
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numerical simulation 2%. These differences may be due to hypotheses and uncertainties in the correlations 
used in the mathematical analysis. Therefore, we observe that the temperature of PV cells with reflectors is 
high due to the increased radiation reflected from the reflectors. 

From the figure 8, we can see the thermal efficiency of PVT collector without Reflectors is near than 
thermal efficiency of PVT collector with reflectors. This result indicates that the reflectors does not have a 
great influence on the performance for this site, because we have different factors that influence the thermal 
performance of the collector such as climatic conditions, incident energy, and outside temperature. The 
geometry of the collector. However, to have a better efficiency, we must maintain a wide temperature 


difference between the inlet and outlet of the heat transfer fluid in addition to raise the thermal conductivity of 
the materials. 
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The electrical efficiency of the PV/T collector with and without reflectors during the day (4 may) for 
the optimal position of the reflectors is presented in Figure 9. In this figure, we can see that the electrical 
efficiency of the collector with reflector is lower than the efficiency without reflectors, due to the increase in 
absorbed solar radiation; this radiation increases the temperature of the cells. 

Therefore, we conclude that the absorbed solar radiation increases the temperature of the cells and 
reduces its efficiency. The electrical yield of the collector with reflector is lower than the yield without 
reflector, which causes an increase in absorbed solar radiation. 
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Figure 9. The effect of reflectors on the electrical efficiency of PVT collector 


4. CONCLUSION 


In this paper, a mathematical model for a solar photovoltaic thermal (PV/T) has been developed. The 
validity of this model was testing by comparing simulation results to experimental ones found in the literature 
where we obtained a good agreement. 

A change of the optimal angle of the reflectors during the day was analytically determined when the 
panel is fixed at B = 35°. Then, we calculated the temperature of the cell and the fluid outlet with and without 
reflectors. Based on the results of the simulation, we found that these temperatures related to the variation of 
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solar radiation during the selected day. On the other hand, the numerical results show that the addition of the 
two reflectors at this site did not greatly improve the thermal and the electrical performances of this collector 
type. 

As a perspective of this work, we aim to improve the performance of hybrid collectors such as using 
nanofluids as cooling fluid because they have high thermophysical properties. 


APPENDIX 
Notation Description Unit 
Ag Collector surface area m? 
Cp Specific heat J/kg K 
m Mass air flow rate kg/s 
M Mass kg 
T Temperature K 
T° Temperature of the previous slice K 
P Power absorbed by glass W 
Nu Nusselt number 
D, Equivalent hydraulic diameter of the useful duct m 
h Convection heat transfer coefficient w/m?.K 
AS Length of a slice along the tube axis m 
At Time of a slice along the tube axis S 
VW, Wind velocity m/s 
Qelec The electrical power produced by the solar cell w/m? 
TSV Local solar time in hours 
N The day number of the day in the year 
i Extraterrestrial intensity on a surface held perpendicular to the sun 
i The solar radiation intensity on a horizontal surface 
Loh The extraterrestrial radiation intensity on a horizontal surface 
Gdir The direct radiation on the collector 
Gdif The sky-diffuse radiation 
Ggr The reflected radiation from the ground 
Gr1 The reflected radiation from bottom reflector which reached the collector 
Surface with titled angle al. 
Gr2 The reflected radiation from the upper reflector which reached the collector 
Surface with titled angle a2. 
Gin The global solar incident radiation 
Gd The global solar diffuse radiation 
Gh The global solar radiation on horizontal surface 
Greek letters 
À Thermal conductivity (W/m K) 
Tg Transmissivity of the glass 
B Fill factor is the ratio of the area covered by the PV cells to the total area of the PV module 
Q The hour angle 
p The latitude of the location 
ô Declination angle 
Subscript 
a Ambient air 
sk Sky 
gr Ground 
W Wind 
I Inlet 
g Cover 
cell Photovoltaic cells 
p Absorber plate 
f Fluid (air Stream) 
pb back plate of absorber 
in Initial or inlet 
O Outlet 
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